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ABSTRACT 

As a radio-loud narrow-line Seyfert 1 galaxy (NLS1) detected by Fermi /LAT in GeV 7-rays, 
1H 0323+342 is a remarkable Active Galactic Nucleus (AGN) showing properties characteristic 
of both NLSls and blazars. Here we present results of simultaneous X-ray and UV/optical 
monitoring observations on 1H 0323+342 taken with the UV/Optical Telescope (UVOT) and 
X-ray Telescope (XRT) onboard the Swift satellite over six years from 2006. Overall, the object 
showed statistically correlated variations in both the UV and X-ray bands on timescales of years 
as well as on timescales of days. A deep Suzaku observation reveals X-ray variability on timescales 
as short as a few tens of thousand seconds, and an X-ray spectrum typical of Seyfert galaxies. The 
broad-band spectral energy distribution, for which the data of UV and X-ray observations taken 
on 2009 July 26-27 were used, can be well modeled with a simple one-zone leptonic jet model plus 
accretion disk/corona emission. The latter is predominantly responsible for the UV/optical and 
X-ray (0.3-10 keV) emission and their observed variations. The correlated UV-X-ray variability 
on the timescale of days is consistent with reprocessing of the X-ray radiation by the accretion 
disk. The shortest timescale and large normalized excess variance of the X-ray variability detected 
with Suzaku suggest a relatively small black hole mass of the order of 1O 7 M 0 , consistent with 
the estimation based on the broad H/3 line in the optical. 

Subject headings: galaxies: active — galaxies: individual (1H 0323+342) — galaxies: perticular — 
galaxies: Seyfert 


1. INTRODUCTION 


As a subclass of Active Galactic Nuclei (AGNs), 
narrow-line Seyfert 1 galaxies (NLSls) are thought 
to have black holes with lower masses and higher 
Eddington r atios compa r ed to normal Seyfert Is 
( Borosoii&Green 1992: Grupe fc Mathur 20041: 


Botte et al. 2004 IXu et al. 20121 see iKomossal 
2008 for a review). While 15%-20% of normal 
AGNs are radio-loud ( Urrv fc Padovanilll995 ). the 
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fraction seems to be much smaller in NLSls. Us¬ 
ing the Sloan Dig ital Sky S urvey Data Release 
3 (SDSS-DR3), IZhou et ah ( 20061 ) built a large 
sample of ~ 2000 NLSls, in which only ~7% were 
detected by the FIRST survey (the Faint Images 
of th e Radio Sky at 20 cm survey; [Beck er et al 


1 19951 ). In a systematic se a rch fo r radio-loud NLS1 
galaxies, iKomossa et al. ( 20061) reported a frac¬ 
tion of 7% radio-loud NLSls, while only 2.5% are 
very radio lo ud ( radio loudness Rsghz > 100). 
Yuan et al.l (120081 ) analyzed 23 genuinely radio- 


loud NLSls with the loudnesf] i?i.4GHz > 100 
and found that some of the objects show obser¬ 
vational properties characteristic of blazars with 


1 The radio loudness f?i.4GHz = /(1.4GHz)//(4400A) in 
lYuan et al.l (120081) is defined as the ratio between the 
flux densities at 1.4 GHz and 4400A in the rest frame, 
instead of commonly used R$g Hz = f(5GHz)/f(4400A) 
iKellermann et, al .1119891) ■ 
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relativistic jets. This finding was confirmed and 
highlighted by the later detections of GeV t 


ray em i ssion with the Fermi/ LA T (jAbdo et al 


2009al lbl: iD’Ammando et 2012 ), as well as the 


rapid variability in the infrared and optical bands 
in three of the extr eme objects, w hic h are also 
detec ted in 7 -rays ( Liu et al. 20101 Jiang et al 


ral structure (IZhou et al . 20071 ) of the host galaxy 
of 1H 0323+342, whereas a ring-like morphology 
was suggested base d on ground-based observations 
( Anton et al.ll2008l) . The asymmetric structure of 
the host galaxy revealed by the analysis of imag¬ 
ing data leads to the suggestion that 1H 0323+342 


brightness temperatures (see, e.g., 

1 7 — n | 

Zhou et al. 

2003 

; Doi et ahl 20061: Zhou et al. 2007t 

Yuan et al. 

2008 

: Giroletti et al. 2011; Doi et al.l 2012) and 


may be associated with a mer ger (jAnton et al 


20081 : Leon Tavares et ahl [20141 ). Based on the 


broad H/3 line width and luminosity, a black hole 


mass o f 1.8 x 10 7 Mq was suggested bv IZhou et al 


the superluminal motion measured from th e radio 
images (see, e.g., ID’ Ammand o et ah 2012i 12013 ) 
also point to the presence of the relativistic jets. 
Given the peculiarity of NLSls among AGNs (i.e. 
relatively small black hole masses, high accretion 
rates, possibly at early phases of the AGN evolu¬ 
tion) , the interesting question arises as to whether 
the jetted NLSls are simply downsized extensions 
of blazars, or whether they form a new population. 

Among the several 7 -ray detected NLSls, 

1H 0323+342 (RA=3 h 24 m 41 s .161, DEC=34°10'45 /, .86, 


(j2007| ). with an u ncertainty of 0.3-0.4 dex as 
commonly inferre d ( Vester gaanl_fc Petersonl[2006l 
Wang et al. l l2009HDong et al.ll 2012 l) . ~ 


At high energies, its 0.3-10 keV spectrum with 
a low signal-to-noise ratio (S/N) obtained by 
Swift/X RT can be fitt ed by a power la w with a 
photon index of T « 2 ( Zhou et aJ1l2007 ). whereas 
the spectrum becomes flatter in the range 20 - 
100 keV, T = 1-55 as obtained by INTEGRAL 
( Panessa et al] 2011 ). Its X-ray emission is highly 
variable (e.g . Zhou et al. 2007 : Abdo et al. 2009bl : 

Paliva et al. 120141) . Th e detection of 1H 032 3+342 


Skrutskie et al. 2006t ) is t he nearest one with a 


Beaalev_et_alJl20Q2 ; also known as 2MASX .103244119+341 (MSB. the Fermi/ LAT ( Abdo et al. 2009bl ) is re¬ 
markable, confirming the presence of a relativistic 
jet. The non-thermal jet emission can account 
for the two bumps of its broad-band spectral en¬ 
ergy distribution (SED), one in the radio-infrared 
band and the other one in the GeV 7 -ray band. 
In additi on, an accretion d isk component is also 
present ( Abdo et al. 2009b ). Interestingly, vari¬ 
ations of the hard X-ray emi s sion above 20 keV 
were found bv iFoschini et al.l ( 20091) : the source 
had a low flux and a soft spectrum as observed 
with INTEGRAL/ IBIS in 2004, and a high flux 
and a hard spectrum as observed with Swift/BAT 
in 2006-2008. This was interpreted as being dom¬ 
inated by emission from the disk/corona during 
the former observation, and by a possible contri¬ 
bution from the jet flaring in the hard X-ray band 
above 20keV during the latter observation. 

1H 0323+342 has been monitored by the Swift 
satellite ( Gehrels et al. 20041) in both the X-ray 
and the UV/optical bands simultaneously over a 
period of more than six years from 2006 to 2013. 
It w as also observed wit h the Suzaku X-ray satel¬ 
lite ( Mitsuda et al.ll2007 ) with a relatively long ex¬ 
posure on 2009 July 26-27. These observations 
provide valuable datasets to study the X-ray and 
UV/optical radiation of this object and their re¬ 
lationship. In this work we analyze the long-term 


redshift of z = 0.0629 (IZhou et all 120071 ). allow¬ 
ing detailed observational studies. It was first 
discovered as an X-ray source in the HEAO-1 
X-ray survey ( Wood et al. Il984 ) and was spec¬ 
troscopically identified with a Seyfert 1 galaxy 
in the optical w ith strong Fe 11 lines and weak 
forbidden lines ( Remillard et al. i fl993l) . Its hy¬ 
brid nature sharing properties of both NLSls 
(e.g., a small broad line width of 1600kms _1 , 
strong Fe 11 emission lin es) an d blazars was first 
reported by Zhou et ahl (2007), making it a pro¬ 
totype of this newly recognized class of AGNs. 
The radio source has a flat radi o spectrum with 


an index ay = +0.1 (S„ oc u ar ) Neumann et al 


1994 iDoi et a l.||2012|). and a h igh radio loudness 
R^art?, = 246 (Doi et ahllioiit ) or i?i. 4 GHz = 318 
dFoschinil 120111) . The VLA 1.4 GHz map re- 
vealed a core with a two-sided structure of ~15kpc 
(lAnton et alJl2008l) . The brightness temperatures 
Tg estimated from both the radio images and flux 
variations are larger than 5 x 10 10 K, exceeding 
the upper limit expected from energy equiparti- 
tion ( Waiima et al. [2014 ). All these observations 


suggest the presence of a relativistic jet in this 
object. The optical image taken with the Hubble 
Space Telescope (HST) revealed a one-armed spi- 
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X-ray and UV/optical data of 1H 0323+342 to 
study both the temporal and spectral properties 
and their variations in these wavebands. Of par¬ 
ticular interest, we find statistical correlations of 
the emission between the X-ray and UV bands on 
both long timescales of years and relatively short 
timescales of days. In Section [2J the observations 
and data reduction are described. The temporal 
and spectral analysis are presented in Section [3] 
and Section U respectively, followed by the broad¬ 
band SED. A summary of the results, and their 
implications are presented and discussed in Sec¬ 
tion [5] Errors on parameters of spectral modeling 
are quoted at the 90% confidence level unless men¬ 
tioned otherwise. Throughout this paper a cos¬ 
mology is assumed with Hq = 70kms -1 Mpc -1 , 
U A = 0.73 and U M = 0.27. 


2. OBSERVATIONS AND DATA RE¬ 
DUCTION 


2.1. The Swift Data 

The Swift satellite has monitored 1H 0323+342 
at 84 occasions from 2006 July to 2013 October, 
with d ata from the X-ray Telescope (XRT, 0.3- 
lOkeV lBurrows et al . 20051 ) and/or t he UV/Optic al 
Teles cope (UVOT, 1700A-6000A, iRoming et al 


20051) . Some of the observations were performed 
at relatively high cadence with timescales of ~1 
day (as shaded in Figure [l]). We have retrieved 
data from the public S wift archives, includin g ob- 
servations publishe d in iFoschini et al.l ( 20091) and 
Paliva et al.1 ( 20141) . All data are reduced using 
the Swift FT00LS in HEAS0FT V6.12. 

The XR T observations in the photon-counting 
mode ('PC’. Illill et al. 2004 ) are used and the data 
are reduced using the task xrtpipeline VO. 12.6. 
X-ray events are selected with grades 0-12 and ex¬ 
tracted using xselect V2.4. To check whether 
the data suffered from the piled-up effects, 0.3- 
10 keV X-ray spectra are first extracted from a cir¬ 
cle of 47" radius (enclosing 90% of the point spread 
function (PSF) at 1.5 keV) centered on the source 
on the observed X-ray images. Some spectra have 
count rates higher than 0.5 counts s - /, indicating 
that they may be affected by pile-upj. To elimi¬ 
nate the effects of pile-up, the central region of the 
source images on the detector are excluded. The 


2 http: / / www. swift. ac. uk/analy sis/xrt/pileup. php 


size of the exclusion r egion is determined by fit¬ 
ting a King function ( Moretti et al.l[2005l) to the 
source image of the data with the highest count 
rate s, and an inner e xclusion radius of 10" (see 
also Abdo et al. is found. Thus a source ex¬ 

traction region avoiding possible pile-up pixels is 
defined to be an annulus with an inner radius of 
10" and an outer radius of 70". For constructing 
light curves, source counts are extracted from such 
an annulus centered at the source position for all 
the datrd; while for extracting source spectra, it is 
only applied to those data suffering from pile-up. 
The background region is chosen to be an annulus 
of the same area centered on the source with an 
inner radius of 75" and an outer radius of 102". 
The spectra are re-binned using grppha V3.0.1 to 
have at least 20 counts per bin. The ancillary re¬ 
sponse function files are created by xrtmkarf and 
corrected for vignetting, hot columns and bad pix¬ 
els using the exposure maps. We use the relevant 
respond matrix given in the output of xrtmkarf. 

During most of the Swift observations, the ob¬ 
ject was also observed with UVOT in all the six 
passbands in the optical (v, b , u bands) and near 
UV (ml, m2, w2 bands). Images of the same 
band during each observation are stacked after as¬ 
pect correction. The source region is chosen to 
be the standard aperture of a circle with a 5" ra¬ 
dius, and the background is estimated from a cir¬ 
cle with a 15" radius near the source region with- 
out contaminatio n from other sources (see also 
Abdo et al. 2010l) . The photon counts are first 


converted into fluxes and then magnitudes in each 
band using the conversion factors with the task 
uvotsource in th e Swift data analysis package 
(IPoolc ct al.ll2008r) , which was obtained by assum¬ 
ing a GRB-like power-law spectrum with dust ex¬ 
tinction. However, the b — v colors derived are out 
of the range valid for this conversion factor, in¬ 
dicating a different spectral shape from the one 
assumed as default. We thus calculate a new 
conversion factor (as well as the effective wave- 


3 We have examined the possible effect of the uncertainty 
of the source positio n (which is better than 3" (~1 pixel, 
iBurrows et all I2005T) 1 on the extracted source counts by 
offsetting the center of the source extraction region by one 
pixel in various directions from the nominal source center. 
It is found that the resulting differences in the extracted 
source counts are at a level of only a few percent at most, 
which are negligible compared with the flux variations of 
the source (see below). 
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length and Galactic extinction) for each of the 
passbands applicable to our target using its actual 
spectral shape. This is done by folding the spec¬ 
trum with the effective area of each of the filters 
of UVOT following Raiteri et al. (2010). The ob¬ 
served spectrum is determined by fitting the aver¬ 
age UV/optical photometric data of 1H 0323+342 
with a log-parabolic model. The Galactic extinc¬ 
tion curve from Fitzpartick (1999) is employed 
with E(B — V) = 0.2095 for this object taken from 
Schlegel et al. (1998). The calibration results are 
shown in Table [T] Then the de-reddened fluxes of 
1H 0323+342 in the UVOT bands are calculated 
using the new conversion factors derived above. 


2.2. The Suzaku Data 


The X-ray Imaging Spe ctrometer (XIS, 0.2- 


ng 

12keV. iKovama et al.l 120071 ) onboard the Suzaku 
satellite observed 1H 0323+342 on 2009 July 
26-27 (ObsID: 704034010, PI: Hayashida, M.) 
for about two days with a total net exposure 
of ~84 ks. Available data are provided by the 
front-illuminated (FI) CCDs XIS0, XIS3, and the 
back-illuminated (BI) CCD XIS1 (XIS2 was lost 
in 2006). We use the Suzaku FT00LS version 19 
and associated CALDB released on 2013 March 5. 
The data are first re-calibrated and screened us¬ 
ing aepipeline. Generally an extraction circle of 
260" radius encircling 99% of a point source flux 
is recommended. In order to keep enough region 
for the background extraction as well as to com¬ 
pensate the lost photons in the wings of the PSF, 
the source light curves and spectra are extracted 
from a source-centered circle of 210" radius, and 
the background light curves and spectra are ex¬ 
tracted from a source-free circle of the same radius, 
avoiding the chip corners containing the calibra¬ 
tion sources. To increase the S/N, we combine 
light curves of all the CCDs, and also light curves 
of observations in both the 3x3 and 5x5 modes, 
into a single light curve with bins of 256 s. For 
each CCD exposure, the data are examined and 


found to be free of pile-up effects (lYamada et al 


2012h by checking the PSF of the X-ray images 
and the count rates at the image peak. The Re¬ 
distribution Matrix Files (RMF) and Ancillary 
Response Files (ARF) are generated by using the 
tasks xisrmfgen and xissimarfgen, respectively. 
Then the spectra from XIS0 and XIS3 are com¬ 
bined by using addascaspec to increase the S/N, 


whereas the spectrum from XIS1 is considered 
separately since the BI CCD has a distinctly dif¬ 
ferent response from the FI CCDs. The spectra 
are grouped to have a minimum counts of 200 per 
bin. Channels in the 1.7-2.3keV range are ignored 
because of the possible uncertain calibration of an 
instrumental Si K edge at ^1.84keV, a I(a fluo¬ 
rescence line at ~1.74keV and an Au M edge at 
~2.2ke’V0. The spectr al analysis is p erformed by 
using XSPEC VI2.7.1 ( Arnaud 19961) and the y 2 
minimization technique is applied. 

We also analyze the data from the PIN silicon 
diodes and GSO phoswich counte rs of the Hard X- 
ray D etector (HXD, 10-600keV. iTakahashi et~ah 
20071 ). The source was not significantly de¬ 


tected above 20keV, however, because of the high 
background. The spectra from PIN and GSO 
are extracted automatically by using the tasks 
hxdpinxbpi and hxdgsoxbpi after recalibration 
and screening. Then we re-bin the spectrum and 
exclude those with negative net counts. The up¬ 
per limit on the hard X-ray fluxes is estimated in 
the following way. The observed number of counts 
Ct in each bin is expected to follow the poisson 
statistics. So the expected number of counts is dis¬ 
tributed as <5(0.9, C t ) at the (one-tail) probability 
of 90%. The upper limit on the source counts 
at the 90% confidence level can be estimated as 
< 5(0-9. Ct) — Ch, wh ere Cb is the background counts 


( Yuan ct al. 20141 ). The counts are converted to 


fluxes using Xspec by assuming a power-law spec¬ 
trum with r = 1.9, similar to that at lower energy 
(see Section Phil) . 


3. VARIABILITY ANALYSIS 
3.1. Long-Term Variability 

The multi-wavelength variability of 1H 0323+342 
is investigated by using the X-ray count rates 
and UV/optical fluxes derived from the Swift and 
Suzaku observations between 2006 and 2013. The 
long-term light curves in various bands, of which 
each data point represents one observation, are 
shown in Figure [T] We test for the presence of 
significant variability at the UV/optical by us¬ 
ing the x 2 test against the null hypothesis of no 
variation. The source is found to be significantly 
variable in essentially all bands with p-values of 

4 http: //heasarc.gsfc.nasa.gov/docs/Suzaku/analysis/abc/ 
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P < 10~ 8 . A trend is apparent that the UV and 
X-ray variability follow each other on timescales 
as long as years, and it becomes less significant in 
the optical bands. To verify this correlated varia¬ 
tions, the UV fluxes, which are obtained from the 
count rates as described in Section [2711 are plot¬ 
ted against the X-ray count rates for the whole 
dataset in Figure[2] Albeit with large scatter, the 
trend of correlated variations between the X-ray 
and the UV bands is evident. This is confirmed 
by the n on-parametric Sp earman Rank correla¬ 
tion test ( Press et aI1ll992ri , which gives extremely 
small values of the null probabilities (p- values are 
P = 2.8 x 1CT 6 , P = 1.3 x 10" 7 and P < 1CT 8 for 
the w 1, m2 and w2 band respectively, see Tabled 
for the observed correlations arising by chance. 

3.2. Short-Term Variability on Timescales 
of Days 

In this work we focus on the short-term 
variations on timescales of days and less in 
1H 0323+342. For this we make use of data from 
the intensive monitoring campaign at high cadence 
performed by Swift/XRT during 2010 October- 
November (shaded region in Figure [T]). Moreover, 
the Suzaku X-ray data taken with a long exposure 
on 2009 July 26-27 are also used. 

For the UV/optical light curves during 2010 
October-November, we first test the significance 
of variations using the y 2 test. The source is con¬ 
sidered to be variable if the chance probability 
P < 0.1%. We find that the UV/optical vari¬ 
ability of 1H 0323+342 is wavelength dependent. 
While no variability is detected in the 6, u, and ul¬ 
traviolet wl band, the object is significantly vari¬ 
able in the v band (with P = 1.2 x 10 -4 ) and 
in the two ultraviolet bands m2 and w2 (with 
P = 3.8 x 10” 4 and P < 10 -8 for the m2 
and w2 bands respectively). Two factors may 
account for this difference. First, it has been 
known that in AGNs the UV/optical variability 
ampl itude increases w i th decreasing wave l ength 
e.g. 


extends out to at least 8 " due to its proxim¬ 
ity and contributes about half of the total opti ¬ 


cal light of the whole galaxy (jZhou et al. 2007]). 


Cutri et al. 1985 : Vanden Berk et al ] 203 


Ai et al.l ' 2013). Second, there is an increasing 


amount of the host galaxy starlight with increas¬ 
ing wavelength towards the optical, which dilutes 
the variations of the central AGN. This is par¬ 
ticularly true for 1H 0323+342 for its prominent 
host galaxy as revealed by the optical HST im¬ 
age taken with the F702W filter (6919A), which 


Such a size is larger than the extraction radius of 
5" used in the photometry analysis above (Sec¬ 
tion m- We check the HST snapshot images of 
1H 0323+342. Unfortunately, we find that the 
images are saturated at the central point source 
(AGN) due to overexposure of the imaging ob¬ 
servations, which hampers any reliable imaging 
analysis of the AGN-host decomposition. How¬ 
ever, we expect that the contamination from the 
host galaxy light in the UV bands is rather low. 
This is because there is evidence that the cen¬ 
tral region of the host ga laxy is dominated by a 
bulge ( Leon Tavares et ai1l2014 ). and consistently, 
our new Palomar 5-m spectroscopic observations 
in the near-IR on this object reveal an old stel¬ 
lar population (Zhou et al. in preparation). We 
thus conclude that the UV fluxes are essentially 
dominated by the AGN flux and use them in the 
variability and SED analysis below. A close-up 
of the UV and X-ray light curves in the time in¬ 
terval 2010 October-November is shown in Fig¬ 
ure [3] For the UV, only the w2 band data are 
displayed, because the sampling with the other 
two filters is rather sparse. Significant variations 
in the UV and particularly the X-ray bands are 
evident. Some interesting features can be seen. 
One may notice that a rapid drop and recovery 
of ~0.4 mag within a single day at the epoch of 
day 30 is observed in the w2 light curve. Simi¬ 
larly rapid changes, of lower amplitudes, occur at 
other time intervals. We have carefully checked 
the images of those datasets for the aspect cor¬ 
rection and other photometric issues, and confirm 
that the variations are real. 

In order to make direct comparison of the vari¬ 
ability between the UV and X-ray bands, we calcu¬ 
late the frac tional variability amp litudes F var for 
these bands (jVaughan et al.l1200311 . 


(f) ' 


(1) 


where a 2 is the variance of the light curve, (e 2 ) 
the mean of the squared measurement errors, and 
(/) the mean flux. The l-cr uncertainty of + var 
due to measurement errors is estimated following 
I Vaughan et al .1 ( 2003 1. The results are given in 
Table [2j showing a much larger fractional vari- 
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ability amplitude in the X-ray than in the UV 
band. Another often used estimator of the vari- 
ability amplitude is the normalized excess v ariance 
( Nandra et al. 1997 : Vaughan et al.ll2003h . which 
is the square of F var and conveys exactly the same 
information as F var . We will talk about this esti¬ 
mation in Section [5.41 


The X-ray variability on even shorter timescales 
is investigated using the Suzaku observations. Fig¬ 
ure Q] shows the Suzaku/XIS light curves in the 
total 0.2-12keV band (upper panel), and in the 
soft (0.2-2 keV, middle panel) and hard band (2- 
12keV, lower panel). The light curves are modu¬ 
lated by the short orbital period for Suzaku 's low- 
altitude orbit (gaps between exposures). It can be 
seen that during the Suzaku observation the source 
varied dramatically by a factor of about two on a 
timescale as short as ~20 ks. The variations seem 
to be independent of the energy bands during the 
Suzaku observation since the same variability is 
seen in both the soft and hard band. 


3.3. Correlated Short-Term UV and X-ray 
Variations 

As indicated in Figure [3l the two light curves 
(w 2 and X-ray) appear to show somewhat coor¬ 
dinated variation patterns. A significant drop fol¬ 
lowed by a rise up within a few days occurred both 
in the w 2 and X-ray bands around day 10 and this 
behavior appears to be simultaneous. However, 
the degree of correlated variations is much lower 
in the second half of the light curve. From day 15 
to day 22, the X-rays varied strongly by nearly a 
factor of two whereas the w2 light curve remained 
stable. After day 25, both X-rays and w2 varied, 
but the correlation between them seems not to be 
so obvious. There seems to be another drop in w2 
from day 23 but its presence is not clear because 
of a large gap between observations. 

Since the data samplings in the two bands are 
nearly paired up during 2010 October-November, 
we can examine any possible statistical correla¬ 
tions by using a correlation test. Figure [5] shows 
the X-ray count rates versus the w2 fluxes for each 
pair of the observations. The results of the Spear¬ 
man correlation test are given in Table [2] As can 
be seen, there is a statistically significant correla¬ 
tion between the two bands, with a null probabil¬ 
ity 1 x 10 -5 . Using the UV fluxes in the other two 
bands (u>l, m2) also yields statistically significant 


correlations between the UV and X-rays (Tabled, 
although the significance is relatively lower due to 
the more sparse data points. 


Given the overall statistical correlations of the 
X-ray and w2 light curves, we search for any possi¬ 
ble time lag between them. The cross-correlation 
function (CCF) method is used, which is com¬ 
monly employed to search for time lags between 
two sets of time ser i es da t a in AGN studie s 
( Hufnagel fe Bregmanl I1992L Kasni et all l2000h . 
To handle the unevenly sampled light curves, 
two complementary methods are adopted here, 
tire i nterpolated cross-correlat ion fun ction ( ICCF ; 
Gaskell fc Sparkdll98fil : iGaskell fc Peterson! 1987 1 


with interpolation bins of 1 day and the z- 
transformed disc rete correlation function (ZDCF; 
Alexander 1997ll . As indicated in Figure [6l both 
methods give generally consistent results and re¬ 
veal a broad signal peaking at r ~ 0 day (pos¬ 
itive time lag means the X-rays are leading the 
UV). We fit the ZDCF in the —6 < t < 9 range 
with a Gaussian profile and find it peaking at 
r ss 0.3 day. A maximu m likelihood method in¬ 
troduced by [Alexander! (20131) gives a peak at 
r = O.Gtfl day with 1 -a uncertainty. These 
results may indicate a correlation with a possi¬ 
ble time lag around zero days between the two 
bands with X-ray tentatively leading. However, 
the uncertainty is large and the lag is statistically 
consistent with zero, given the broad peak of the 
CCF. 


To test the significance of the cross-correlation 
function, we simulate X-ray light curves that are 
uncorrelated with the observed one and then cal¬ 
culate the ICCFs between simulated X-ray light 
curves and the observed UV light curve. The de¬ 
tailed procedure of the simulations is described in 
Appendix [A] The 68%, 95% and 99% extremes of 
the ICCFs distribution out of the two uncorrelated 
light curves are indicated by dashed, dash-dotted 
and dotted lines in Figure 0 As shown in the 
figure, the peak around zero lag in the real data 
reaches a significance level just slightly below 95%, 
and thus can only be considered to be marginal. 
The lower significance level found for the corre¬ 
lation than that from the direct correlation anal¬ 
ysis using the Spearman’s Rank correlation test 
(Table [2]) is mostly due to the random fluctua¬ 
tions introduced from the interpolation of the two 
light curves by assuming that the light curves vary 
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smoo thly fe.g. lWhite &; Petersodl 1994 1 Alexander 
ll997 ). We note that there appears to be a second 
peak around r ~ —10 days. However, given the 
even lower significance and the shortened segment 
of data involved in the correlation (~ 25 days), 
we do not consider this peak to be a real corre¬ 
lation, rather a coincidence of the second UV dip 
(~23 day, see Figure EJ) with the first X-ray dip. 
In other parts of the whole long-term light curves, 
the number of paired data points with good sam¬ 
pling cadence is too small to carry out meaningful 
cross-correlation analysis. 

4. SPECTRAL ANALYSIS 
4.1. X-ray Spectral Analysis 

The X-ray spectra observed with Suzaku/XIS 
are analyzed for their much higher S/N compared 
to the Swift/XRT spectra. The combined FI spec¬ 
trum in the band 0.5-10 keV is used. We first con¬ 
sider the total, time-averaged spectrum, as shown 
in Figure[ 8 ] All the best-fit parameters are listed 
in Table [3] An absorbed single power-law model 
with the absorption fixed at the Galactic value 
(A^ al = 1.27 x 10 21 cm -2 ; iKAlherla et, a~i1l2005l) 
gives an unacceptable fit (y 2 /d.o.f.=923/362), 
with a systematic excess of flux below 1 keV in 
the residuals (Figure [ 8 ] b). A black body (zbbody, 
with a redshift fixed at z = 0.0629) is then added 
to account for this soft excess, which improves 
the fit significantly (y 2 /d.o.f.=354/360), giving a 
black body temperature kT = 0.15 ± 0.01 keV and 
a photon index T = 1.87 ± 0.02. The temper¬ 
ature is within the previously observed range ~ 
0.1-0.2keV of the soft X-ray excess in t he X- 
ray spectra o f AG Ns (i Gierliriski fe Done] 2004t 
Porquct et ah 20041 Crummy et al.l2006 : Ai et al 


rqu 

201111 . The photon index is somewhat flat for 
NLSls but not unseen, because there is a large 
scatter in th e distribution o f spectral ind i ces of 
NLSls (e.g 


Lcighlv 1999bi : Grupe et al.l l20ld : 


Ai et al.l l201lh. We also consider a double power- 


law model, which gives a good fit. Alternatively, 
to account for the soft X-ray emission, we fit a 
partial covering model with io nized absorption 
(zxipcf in Xspec, Beeves et al. 2008t h which re¬ 
sults in a somewhat worse y 2 and a marginally 
acceptable fit with a lowly ionized absorber of 
Ah = 6 x 10 22 cm -2 and a moderate covering frac¬ 
tion of / = 58%. We also fit a Comptonization 


model (comptt in Xspec. lTitarchuldll994r) and fix 
the seed photon energy at the innermost temper¬ 
ature of the standard accretion disk based on the 
black hole mass and accretion rate of this object 
(see below). The model results in an improved fit 
over the black body model, and the fitted elec¬ 
tron temperature fcT p i asma = 0.28 keV and optical 
depth r = 13.4 are within the previou sly observed 
range in similar ob jects (e.g. Gierliriski fc Donel 


20041; lAi et ai]|201lh . 


In addition to the soft X-ray excess, there is 
a systematic excess in the residuals above 6 keV 
(Figure[ 8 t,d). We have also used a primary power- 
law continuum plus a reflection com ponent. The 
i onize d disk reflection model from lRoss fc Fabian 
(|2005l . reflionx in Xspec) is employed. In or¬ 


der to account for the relativistic blurring caused 
by strong gravitation near t he black hole, we use 


relconv ( Dauser et al. 2010h as a relativistic con¬ 


volution kernel to convolve with the disk reflec¬ 
tion component. The free parameters are the in¬ 
dex of the primary power-law component T, the 
black hole spin a, the abundance of iron rela¬ 
tive to the solar value Ap e and the ionization 
pa rameter £r. A simi lar approach was adopted 
bv lWalton et al. ( 2013 ) to investigate the Suzaku 
observations of a sample of AGNs including 1H 
0323+342. They found that the inclination ob¬ 
tained was unphysically high for an unobscured ac¬ 
tive galaxy (i = 82° ± 3°), so the inclination angle 
was fixed at 45° in the fitting. Here we fixed the in¬ 
clination at 10 °, considering the small inclination 
angle of the jet inferred from its blazar properties 
and detected 7 -ray emission. The inner radius of 
the accretion disk is fixed at the innermost stable 
circular orbit and the outer radius fixed at 400r g , 
where r g is the gravitational radius. The radial 
emissivity profile of the disk is assumed to be a 
power law e oc r~ q with q = 3. The model gives 
a very good fit (y 2 /d.o.f.=351/355). The spin pa¬ 
rameter, though poorly constrained, is found to be 
< 0.13 at the 90% c onfidence level . This result is 
different from that of I Walton et al.l ( 20131) possibly 
due to the very different inclination angle adopted 
here, but consistent with the fin ding of generall y 
low black hole spins for NLSls ( Liu et akl 20 l|). 
However, a reliable estimate on the black hole spin 
from fitting the reflection component can only be 
achieved by extending the observed X-ray spec¬ 
trum to higher energies up to 30-40 keV. 
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Fitting the spectrum from the single BI CCD 
(XIS1) gives essentially consistent results with 
those for the FI CCD spectra. Since adding the 
single BI spectrum results in only limited improve¬ 
ment in the S/N, we do not incorporate it in the 
analysis, considering possible systematics in the 
cross-calibration. 

The light curve measured with Suzaku reveals 
strong variations up to a factor of two on short 
timescales. To investigate possible spectral evolu¬ 
tion with changes in the flux, we extract the FI 
spectra at the epochs of the highest and the low¬ 
est count rates, respectively, in the Suzaku light 
curve in Figured] (labeled by bars). Each spec¬ 
trum has a minimum of 50 counts per bin. The 
high- and low-flux spectra show very similar pro¬ 
files, as indicated in Figured] They are then fitted 
with the above power law plus black body model, 
respectively. Although the high-state spectrum 
has a slightly steeper best-fit index value than the 
low-state, the fitted parameters of the two spectra 
are consistent within the mutual uncertainties, as 
listed in Table dl We also calculate the flux ra¬ 
tios of the black body to the power-law compo¬ 
nent within 0.3-10 keV for the two spectra, and 
find them to be the same. Then we follow a dif¬ 
ferent approach. We extract the high-state and 
low-state spectrum by separating them above and 
below the count rate threshold of 2.4 counts s -1 , 
and fit them with the above model. Both spectra 
have photon indices of T = 1.87 and black body 
temperature of kT = 0.15 keV. The flux ratios of 
the black body to power-law component for high 
and low state are still the same. These results are 
consistent with the energy-independent variabil¬ 
ity as seen in Figure 0] We thus conclude that, 
on timescales as short as tens of thousand seconds 
and with the variations within a factor of two or 
so, there is little or no noticeable change in the 
spectral shape. 

To investigate possible spectral variability on 
timescales of days, the Swift/XKT spectra in the 
0.3-10 keV band are extracted from each of the ob¬ 
servations taken during 2010 October-November, 
as described in Section 12.11 Although there is 
a soft excess in its X-ray spectrum, the spectral 
S/N are generally low, and the purpose here is 
to find the variation of the spectral shape possibly 
caused by, e.g., changes of the ratio of the soft and 
the hard components. So we fit these spectra us¬ 


ing a single power-law model with the absorption 
column density fixed at the Galactic value. The 
best-fit photon indices are plotted in Figure [10] 
No spectral variability is detected (P = 0.2 using 
the x 2 test). 


4.2. Broad Band SED 

Given its hybrid nature, sharing properties of 
both NLSls and blazars, it would be revealing 
to stu d y th e broad band SED for 1H 0323+342. 
lAbdo ct~ahl ( 2009bh reported an SED using the 
averaged X-ray and UV/optical data from the 
Swift observations between 2006 July and 2008 
November. It showed a low synchrotron-peak fre¬ 
quency at about 10 13 Hz and a strong disk/corona 
component. Given the relatively large variability 
of the source, it would be ideal to construct an 
SED based on simultaneous observations. Here 
we make use of the X-ray data from the Suzaku 
observation taken on 2009 July 26-27 and the 
UV/optical data taken from the Swift obser¬ 
vation on 2009 July 27, which can be consid¬ 
ered as quasi-simultaneous. However, the data 
in the other wavebands were not simultaneously 
taken. The GHz radio data are taken from the 
NASA/IPAC Extragalactic Database^ and listed 
in Table[5] Data in the infrared, hard X-ray and 7 - 
ray band s are retrieved f rom the A11WISE Source 
Catalog ( Cutri fc ct alJ 2012 1. t he Swift /BA T 70- 


Month Catalog ( Baumgartner ct al. 20131 ) and 


the Fermi /LAT 2-year Point Source Catalog 
( Nolan et al l lioii . respectively. The broad band 
SED is shown in Figure HU 

The broad-band SED is fitted with a sim¬ 
ple one-zone leptonic jet model, which con¬ 
sists of synchrotron, synchrotron self-Compton 
(SSC) and external-Compton (EC) processes (e.g . 
Ghisellini fc Tavecchiol 120091 : IZhang et alJ 20141 : 


Sun et all 20141) . plus emission from an accretion 


disk/corona. The % 2 minimization techni q ue is 


appli ed to perform the fitting (see IZhang et al 


I2014L for details). The energy distribution of the 
injected relativistic electrons is assumed to be a 
broken power law in the range of [ 7 min, 7 max] with 
a break energy 7 b. The 7 ma x is usually poorly 
constrained and does not affect the results signif¬ 
icantly, so it is fixed at a large value. The break 
energy 7 b and the energy density parameter Nq 


5 http://ned.ipac.caltech.edu/ 





































can be replaced by the peak frequency of the syn¬ 
chrotron component v syn and the corresponding 
peak lum inos ity i/ SV n l v (see Equation (2) and (5) 
in Zhang et al. 2009h during the fitting, so we set 


r'syn and v syn l u as free parameters instead of 7 b 
and Nq. The two indices pi, P 2 of the broken 
power law are obtained from the spectral slopes 
below and above the s econd bump o f the observe d 


SED as described in Zhang et al. (2012, 2014 1. 


Here the jet radiation region is assumed to be a 
sphere with a size determined by the vari ability 


spnere witli a size determined by tne vari abimty 
of jet emission as adopted by ISun et ah ( 2014j) 
specifically for 1H 0323+342. 

One of the key factors to determine the rel¬ 
ative importance of different sources supplying 
seed photons for inverse Compton is the distance 
of the dissipation region in the jet from the black 
h ole, w hi ch is st ill co ntrov ersial in AGNs (e.g. 
Ghisellini fc Tavccchiol 20091 ). Some of the stud¬ 


ies assume that the process takes place within 
the broad line region (BLR) and the seed pho¬ 
tons are dominated by radiation from the BLR 
(hereafter IC/BLR) because of its much hig her en ¬ 
ergy density than that of the torus (e.g. iBai et al 


120091 : llsler et alj|2013t Izhang et al.ll2014 ), whereas 

in other cases the process is assumed to takes 
place outside the BLR and the seed photons come 
predominantly from a dusty torus as infrared 
emission (IC/IR) since the photon field from the 
BLR will decrease ver y quickly in this case (e.g . 


Tavecchio et al.1 12013 : D’Ammando et all 20121 


20131) . The latter is usually invoked to explain the 


detection of very high energy gamma-rays extend¬ 
ing to the TeV regime in some blazars, which can¬ 
not be accounted for by the IC/ BLR p r ocess due 
to th e Klein-Nishina effect (e.g. iTavecchio et al. 
For 1H 0323+342, no TeV 7 -ray emis- 
95% upper limit of 


20111 ). 


1H 0323+342, 
detected with a 


sion is 

5. 2 x 10 ~ 12 ergs -1 cm ~ 2 at energies above 400 GeV 


(IFalcone et al.H2004+ However, this can not rule 
out the IC of seed photos from a putative torus. 
Here we consider both cases that the dissipation 
region is inside and outside the BLR, respectively. 
In the first case, the dissipation region is inside the 
BLR so that IC/BLR dominates and IC/IR can be 
ignored. The energy density of the BLR in the co¬ 
moving frame, I/ BLR = 1.55 x 10 _ 2 T? et erg cm -3 , 
is obtained from its strong broad Balmer lines 
(~ 14.5 x 10~ 14 ergs -1 cm -2 , Zhou 2007, private 
communication), calculated using the method of 


Sun et al.l (|2014l) . where the Lorentz factor Tj e t 
of the jet is taken as its Doppler boosting factor, 
i.e., Tjet = S, by assuming that its relativistic jet 
is inclined close to the line of sight given its 7 - 
ray detection. Alternatively, we also consider the 
case where the IC is dominated by seed photons 
from a putative dusty torus (IC/IR). The energy 
density U IR = 3 x 10 _ 4 T 2 et erg cm -3 is adopted 
as a conventional v alue, as in previous works (e.g. 
Idearv et al.ll2007h . As shown below, both cases 
lead to satisfactory fits to the non-thermal emis¬ 
sion of the broad-band SED. 

Given the prominent broad emission lines in the 
optical spectrum as in other NLSls, a strong ion¬ 
izing continuum emission, which originates from 
an accretion disk around the black hole, must 
be present in 1H 0323+342. We use a multi- 
tcmp erature black body model of a standard thin 
disk rtDavis fc Laoii l2011r) to describe this compo¬ 


nent. A similar component was also adopted in 
modeling the SED of 1H 0323+342 in the previ¬ 
ous studies (lAbdo et al. 2009b : Paliva et al. 20l4 
ISun et aTI 12014 ). A black hole mass of 1.8 x 
10 6 7 * * * * Mq d etermine d from the broad H/3 luminos¬ 
ity by I Zhou et al.l ( 20071) is employed. The inner 
and outer radii of the disk are fixed at 37? sc h and 
700i? S ch respectively (R sc h is the Schwarzschild ra¬ 
dius) and the inclination is taken to be cosi = 1 . 
The accretion rate is varied to match the UV 
fluxes. It should be noted that the optical fluxes 
can only be considered as upper limits due to pos¬ 
sible contamination from the host galaxy (marked 
as arrows in Figure fill) , while the UV fluxes are 
not (see Section 3.2). For the X-ray (Suzaku) 
data below 10 keV, the best-fit physical model of a 
power law plus Comptonization by warm plasma 
for the soft X-ray excess is adopted (see Section 
4.1). 

We ignore the thermal emission in the infrared 


6 Given the sparse data in the optical to EUV band, we sim¬ 
ply use the multi-temperature black body model of the ge¬ 

ometrically thin, optically thick standard accretion disk to 
fit the SED. It should be noted that this model may not be 
appropriate for 1H 0323+342 because the derived nominal 

value of the Eddington ratio of this object is higher than 

the critical value for the geometrically thin approximation 

to be valid. Thus the fitted Eddington ratio value should 
be treated with caution. Moreover, the large uncertainty 
in the black hole mass based on the single-epoch spectrum 

estimation (~ 0.3 dex) introduces a further uncertainty in 

the derived Eddington ratio. See discussion in Section 5.1. 
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band from a putative dusty torus in considera¬ 
tion of a few arguments. First, in a few simi¬ 
lar radio-loud NLSls, the infrared emission in the 
WISE bands is highly variable on the timescales 
of hours, and is considered to be do minated by 


non-t hermal emission from the jets (IJiang et al 
2012) • Second, based on previous SED model¬ 


ing for several 7 -ray detected radio-loud NLSls 
including 1H 0323+342, the contribution to the in¬ 
frared from an assumed torus is found to be much 
smaller compared with the s ynchrotron radiatio n 
of the jets and is negligible ( Abdo et alJ[2009b ). 
Last, adding a contribution in the infrared from a 
putative torus, which is uncertain on its own any¬ 
way, will have essentially little or no effect on our 
conclusions. 

During the fitting procedure, the free param¬ 
eters .B, 8, v syn , v S ynK and 7 min are constrained 
with the \ 2 minimization technique, while all 
other parameters are fixed at the given values as 
discussed above. The best-fit model parameters 
are listed in Table [ 6 ] for both the two cases, to¬ 
gether with their 1-cr errors. We note that the UV 
fluxes are dominated by the thermal disk emis¬ 
sion and the X-ray band below lOkeV is domi¬ 
nated by the corona emission in both the IC/BLR 
and IC/IR cases, which give satisfactory fits to 
the SED in the infrared, and from the hard X-ray 
(tens of keV) to the GeV 7 -ray band. The rela¬ 
tive contributions of the disk/corona and the jet 
components in the UV/optical and X-ray bands 
are consistent in these two cases. We plot the 
fitting results of the IC/BLR case as a demon¬ 
stration in Figure fill As indicated in the figure, 
only in the hard X-ray band (> lOkeV), the jet 
emission starts to make a noticeable and gradually 
increasing contribution, and eventually dominates 
the observed flux at energies around lOOkeV and 
above. We also note that the infrared data are not 
matched well by the model, except for their over¬ 
all intensities. This may be due to a contribution 
from the host galaxy in the near-infrared. 

Since the UV/optical and X-ray (< lOkeV) 
data are not simultaneous with observations in the 
other bands, the SED may be somewhat differ¬ 
ent at different epo chs given the X-ray variability 
seen in this source ( Foschini et alJl2009|h Here we 
examine to what extent the above modeling will 
differ, due to the source variability, especially con¬ 
cerning the UV/optical and X-ray bands presented 


in this work. We vary the hard X-ray flux in the 
energy band of Swift /BAT (14-195 keV) by a fac¬ 
tor of two in both directions, and re-do the fitting. 
These fitting results are over-plotted in Figure fill 
with grey lines. We find that such changes have 
almost no effect in the UV/optical regime, and the 
UV fluxes are always dominated by thermal emis¬ 
sion from the accretion disk. In the X-ray band 
below lOkeV the contribution from the jet emis¬ 
sion is still negligible compared with the dominant 
corona component. It should be noted that the 
dominance of the thermal disk and corona emis¬ 
sion in the UV and X-rays below lOkeV, respec¬ 
tively, is essentially not affected by varying the rel¬ 
ative contributions of the coronal X-ray luminos¬ 
ity, the black hole mass, and the jet parameters 
within the allowed space that match the model 
SED with the data observed. 


Using our SED model, we can quantify the 
relative contributions from each of the compo¬ 
nents to the observed luminosities in the Swift and 
Suzaku bands. The UV luminosity is integrated 
from 1300A to 330C)A, covering the w avebands of 
the three UVOT ultraviolet filters ( Poole et al„ 
200811 . yielding L uv « 1.9 x 10 44 ergs’ 1 . The 


jet component contributes only 1.3%, Luv,jet ~ 
2.4 x 10 42 ergs _1 . The total 0.3-10keV lumi¬ 
nosity is Lx ss 2.5 x 10 44 ergs _1 , where the jet 
emission from inverse-Compton contributes only 
Lx, jet ~ 1-2 x 10 43 ergs’ 1 , i.e. ~ 5%. The frac¬ 
tion may be even smaller if the soft X-ray excess 
component is taken into account. We thus sug¬ 
gest that, in the context of the simple one-zone 
leptonic model, the jet contribution to the X-ray 
band (<10keV) is most likely negligible compared 
with the disk/corona emission at the epochs of the 
Swift and Suzaku observations considered here. 


5. SUMMARY AND DISCUSSION 

In this work, we analyze the X-ray and UV/optical 
data from the long-term monitoring with Swift of 
the NLSl-blazar hybrid object 1H 0323+342. We 
also analyze the Suzaku observation of this source 
taken on 2009 July 26-27. The main results are 
summarized as below, and their implications are 
discussed in the following subsections, (i) The ob¬ 
ject is variable in both the UV and X-ray bands on 
timescales from days to years, and there is a statis¬ 
tically significant correlation between the UV flux 
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and X-ray count rates on these timescales. On 
timescales as short as a few tens of thousands sec¬ 
onds, X-ray variability by a factor of two is also 
found, though no significant spectral variation 
is detected, (ii) Using observations of intensive 
(roughly daily) monitoring with Swift lasting for 
~35 days, a cross-correlation analysis suggests a 
possible time lag around r ~ 0 day between the 
UV and the X-ray emission with X-rays tenta¬ 
tively leading. However, the present data do not 
allow any firm conclusion, and future monitor¬ 
ing observations are needed to confirm this result, 
(iii) The high S/N X-ray spectrum of 1H 0323+342 
taken with Suzaku/XIS is typical of Seyfert Is, a 
power law with a photon index of T ~ 1.9 and a 
soft excess, the latter contributing to about 11% 
of the flux within 0.3-10 keV. (iv) The broad-band 
SED is constructed with the (quasi-) simultane¬ 
ous UV/optical and X-ray data, as well as the new 
detection of this object in the 100-300 MeV band 
with Fermi/ LAT. The latter is important for pin¬ 
ning down the exact position of the high-energy 
tail of the y-ray bump (thus the energy budget 
of the jet emission), improving upon the loosely 
constrained limit adopted in previous SEDs. (v) 
The SED can be well modeled in the context of 
a simple one-zone leptonic jet model plus the ac¬ 
cretion disk/corona components. The latter are 
suggested to dominate the emission in the UV and 
X-ray (up to lOkeV) bands as observed with Swift 
and Suzaku presented here. 


5.1. A NLSl-Blazar Hybrid and the UV/X- 
ray Emission 

The y-ray NLSls discovered so far are all 
suggested to be low synchrotron-peaked blazar s 
( Abdo et al.l l2009b ; D’Ammando et al. 2012 b 
Our SED fitting reveals that the synchrotron com¬ 
ponent of 1H 0323+342 pea ks at ~ 10 13 Hz, con ¬ 


sistent with this assertion (lAbdo et al.1 l2009bO . 
Considering the strong broad Hfl line with a n 
equivalent width of 58A ( Zhou et al. 2007 ). 
1H 0323+342 is similar to flat-spectrum radio 
quasars (FSRQs), whose synchrotron peaks are 


also at low frequencies in t he infrared (lAbdo et al 


20 ld ; Giommi et al.l[2oT3) . The SED fitting gives 
a bolometric luminosity of Lboi = 1.9x 10 45 ergs -1 
(Tabled here the bolometric luminosity is the sum 
of the disk and corona component). Adopting the 
mass of the black hole 1.8 x 10 7 Mq estimated 


from the H/? line (jZhou et al . 20071 ). the Edding¬ 
ton ratio is then 0.8 based on the SED modeling^. 
It should be noted that Zhou et al. (2007) gave 
Uboi = 1-2 x 10 45 ergs -1 , a factor of 1.6 lower, by 
us ing the bolometric correction factor for quasars 
of lElvis et al.l ( 1994h . By using the UV spectral 
slope cqjv = 1.13 fitted from the three UV fluxes 
and applying the rela tion of cmv and the Edding¬ 
ton ratio for NLSls in drupe et al.l ( 2010h . we find 
a value of 0.51; however, there is a large scatter in 
this relation. Such a high Eddington ratio is char¬ 
acteristic of typical NLSls, and is most likely the 
underlying driver for the strong Fe 11 emission and 
other properties found in 1H 0323+342, as in other 
NLSls. Increasing the black hole mass will result 
in a decreased Eddington ratio, and vice versa. 
However, if the observed extreme properties of 
NLSls are indeed driven by their high Eddington 
ratios as commonly thought, the black hole mass 
of 1H 0323+342 cannot be much higher (e.g., a 
mass ten times higher will lead to an Eddington 
ratio < 0.1). 

The dominance of the thermal disk/corona 
emission in the UV and 0.3-10 keV X-ray bands is 


also consistent with previous results (lAbdo et al 


2009bl) . The jet emission contributes only a few 


percent to the total luminosity in these two bands 
(Section l4.2L though its contribution may increase 
toward higher energ ies in the hard X-ray band 
( Foschini et al] 12009t h It should be noted that, in 
some other y-ray detected NLSls, the X-rays be¬ 
low 10 keV could be dominated by the jet emission, 
at least at some occ asions (e.g., PMN J 0948+0022 
and PKS 2004-447, Abdo et aiJl2009bh . This may 
be caused by the emergence of an inverse-Compton 
component in the X-ray band during the jet flares, 
which seems not to be the case for 1H 0323+342 
at the epochs of the Suzaku and Swift observations 
analyzed here. 


5.2. On the Long-Term X-ray/UV Varia¬ 
tions 

Thanks to the extensively long monitoring ob¬ 
servations carried out by Swift , statistically corre- 


7 We note that such a value should be treated with caution 
since, apart from the large uncertainty in the estimated 
black hole mass (~ 0.3 dex for estimation based on single¬ 
epoch spectra), the actual accretion fl ow may differ from 
the s tandard thin disk model (see, e.g.. lAb ramowic z et all 
1 1988H . which is used to fit the UV data here. 
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lated variability between the UV and X-ray bands 
are found in 1H 0323+342, which are significant 
albeit with large scatter. The long-term flux vari¬ 
ations in the optical/UV in AGNs have been stud¬ 
ied extensively over the past decade s , for individ¬ 
ual objects fe.g. lEdelson et al. 1996 : Mason et al 


20021 : Alston et al. 2013) and for large samples 


(e.g. Cutri et al.1 I l985l: Tre vese fc Vagnettil 2002 : 


20041 Ai et al.1 2013t h It has 
also been found that, in general, NLSls vary 
with systematically smaller amplitudes compared 
to AGNs with larger broad line widths (Ai et 
al. 2013). The long-term UV flux variations in 
1H 0323+342 (Tabled have amplitudes that agree 
well with the median values for the NLS1 sam¬ 
ple in Ai et al. (2013, their Figure7), and thus 
likely have the same origin. If the UV bands are 
dominated by the disk component during the Swift 
monitoring, as suggested by the SED above, the 
long-term UV variability can be well explained as 
variations of the accretion rate, as commo nly be¬ 
lieved to operate in AGNs (e.g. Li &; Caol 2008l h 
The timescale of such variations, in the frame¬ 
work of the standard optically thick disk model 
for an accreting black h ole, is charac terized by the 
viscous timescale r v i S c ( CzernvIfeoOG ). 


Vanden Berk et al 


Tvisc = T th (R/H) 2 = ( rta)-\R/H ) 2 , (2) 

where T t h is the thermal timescale, U the Kepler 
angular velocity, a the viscosity parameter and 
R/H is the ratio of the disk radius to the scale 
height. By assuming a = 0.1, R/H = 10, and 
adopting Mbh = 1-8 x 1O 7 M 0 , we estimate a vis¬ 
cous timescale at R = 10i? sc h of r v i sc = 1 x 10 7 s « 
120 days. The long-term UV/optical variations as 
observed with Swift (Figure [U may be caused by 
such fluctuations in the mass accretion rate. The 
X-ray emission from the hot corona is ultimately 
powered by accretion, though the exact physical 
mechanisms to transfer the accretion energy from 
disk to corona is not understood yet. In any case, 
the processes must be efficient and the timescales 
are likely much shorter than the viscous timescale 
in the disk, suc h as via the p rocess of magnetic re¬ 


connection (Liu et al. 200311. In such a scenario 


I' 


the correlated X-ray and UV/optical variability 
are expected naturally on the timescales of sev¬ 
eral months to years. 


5.3. On the Short-Term Correlated X- 
ray/UV Variations 


Now we consider shorter timescales of about 
several days, on which 1H 0323+342 also varied 
in the X-ray and UV bands (Figure[3] and Ta- 
ble[2]). For the X-ray emitting corona (geomet¬ 
rically thick, R/H ~ 1), the viscous timescale 
i+ sc is comparable to the thermal timescale T th. 
which is much shorter than the viscous timescale 
of the disk. For the above black hole mass, 7 + sc ~ 
r t c h « 10 5 s. Thus the X-ray variability can be in¬ 
terpreted to arise from variations in the corona 
mass density and temperature. The former may 
be related to the fluctuations in the accretion flow 
of the corona gas, and the latter may be related 
to the fluctuations in the energy dissipation (heat¬ 
ing) process of the corona. It is thus expected for 
the X-rays to vary on timescales of days or longer. 

The UV variability on these timescales can¬ 
not be explained by fluctuations of the mass 
accretion rate, since these timescales are much 
shorter than the viscous timescale of the disk. 
Such variability is commonly seen in AGNs from 
the moni toring observations for individual ob¬ 
jects (e.g. Edelson et al. I ll996t iMason et al.11200'3 : 


Alston et al. 2013 1 as w ell as for large samples 


(e.g. Collier fc Petersoi] 12001 1. On average, on 
timescales of days, the UV/optical variability am¬ 
plitudes are small (typically several percent) for 
Seyfert galaxies; they are even smaller for NLSls 
based on their ensemble structure functions (Ai et 
al. 2013, their Figure 10), which show a rise as a 
power-law function with a slope of 0.3-0.4. The 
UV variability amplitudes found in 1H 0323+342 
during the 35-day intensive monitoring (Figure 3 
and Tabled are roughly comparable to those of 
the Seyfert 1 population. The physical mechanism 
governing the short-term UV/optical variations is 
not well understood yet, however. 


One of the most interesting results of this work 
is the evidence for statistically correlated UV- 
X-ray variability on timescales of days. This 
is observed for the first time among 7 -ray de¬ 
tected NLSls. Similar correlated UV/optical- 
X-ray variations on similar timescales have been 
previously reported in a limited number of AGNs, 
mostly radio-quiet, with intensive simultane¬ 


ous monitoring observations, e.g., N GC 4051 
( Breedt et al. 201fll : Alston et al. 2013lh MR2251- 
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Arevalo et al. 2008), Mrk79 ( 

Breedt et al. 

20091). PC 1211+143 (Baehevetal. 

2009). The 


most promising explanation invokes reprocessing 
of the primary X-ray emission by an irradiated 
‘cold’ accretion disk, whose local temperature 
may, in turn, be altered by the reprocessed en¬ 
ergy dissipated into the disk. Changes in the ir¬ 
radiating X-ray emission would cause subsequent 
changes of the local emissivity in the disk, leading 
to the cor responding variatio ns in the UV/optical 


emission (Colli n ct al.| 120031 ). The reprocessing 


model has been used to explain the UV/optical 
variations on the timescales of days and their 


AGNs (e.s. Clavel et al. 1992: 

Nandra et al. 

1998 

200C: Shemmer et al.l 2001: 

Mason et al. 

2002 

Arevalo et al. 20081 Breedt et al. 2009, 

201 C 

Alston et al.| 2013). 


Based on the argument of energetics, the 
changes in the UV/optical luminosity should al¬ 
ways be smaller than (or at most comparable to) 
the changes in the X-ray luminosity, which is the 
driving component. In the case of 1H 0323+342, 
during the 35-day monitoring in 2010 October- 
November, the amplitude of the luminosity varia¬ 
tions in X-ray band is A Lx = 6.1 x 10 43 ergs _1 , 
compared to ALuv < 1-9 x 10 43 ergs -1 in the UV 
band. The fractional variability amplitude in the 
X-ray band, F var ,x = 29%, is also larger than that 
in the UV, F V ar,uv < 10%. 


When a variation in the driving X-ray emission 
occurs, a time lag is expected for the UV/optical 
emission to respond, which is roughly the readjust¬ 
ment timescale (dynamical timescale) of the sur¬ 
face layer of the disk, heated by the irradiating X- 
ray emission. Such a timescale is generally short. 
As discussed in Collin et al. ( 20031 ). assuming that 
the Thomson thickness is unity and the readjust¬ 
ment takes place at the sound speed, this timescale 
for an AGN with 10 s Mq black hole mass is: 


t dyn ~ lO 5 ^- 172 ^ 1 s, (3) 


where Tq is the disk temperature in units of 10 6 K 
and ni 2 the electron number density in units of 
10 12 cm -3 . For 1H 0323+342 with a smaller black 
hole mass, this timescale could be < 10 5 s. The 
time lag of around zero days with a possible lead¬ 
ing by the X-rays, as indicated by the cross¬ 
correlation analysis in Section 13.31 is consistent 


with this estimat(jf|. More intensive monitoring 
with even higher cadence of data sampling than 
the current observations is needed to confirm this 
result, however. 

Alternatively, in light of the fact that the ob¬ 
served time lag is consistent with zero, the corre¬ 
lated X-ray and UV emission may also be caused 
by the following effect. As one can see from Fig¬ 
ure 4, the variation in the full 0.2-10 keV band 
is dominated by that in the soft X-rays in 0.2- 
2keV, where the soft X-ray excess dominates. It 
has been shown in Section 4.1 that the soft X- 
ray excess can well be modeled as the high-energy 
tail of Comptonized emission of seed photons from 
the underlying accretion disk by warm plasma of 
0.28 keV with an optical depth of r ~13.4 (i.e. 
the mean number of scatterings = r 2 ~ 180). If 
this is indeed the case, the observed UV emission 
may also be contributed by Comptonized photons, 
which arise either from seed photons of low ener¬ 
gies and/or from experiencing only a small number 
of scatterings. Fluctuations in the warm plasma 
may give rise to synchronized variations in both 
the (soft) X-ray and UV bands naturally with no 
time lag. 

In reality, the interplay between the X-ray and 
UV/optical variations must be complex rather 
than a simple one-to-one follow-up, considering 
the feedback of the UV/optical radiation to the 
X-rays as seed photons for Compton scattering, 
as well as other possible factors. For example, the 
observed variations of the X-rays may be caused 
by small changes in the ionization state or the 
covering factor of a (dust-free) ionized absorber, 
whereas the UV/optical emission would be left un¬ 
affected. This may explain the mis-match between 
the X-rays and UV seen in the latter half of the 
light curves in Figure [3] 

Finally, based on our SED modeling, the jet 
contributes only < 5% to the UV and 0.3-10keV 
X-ray band. Therefore, though the jet contribu¬ 
tions cannot be ruled out completely in the hard 
X-ray band, it is unlikely that the observed varia¬ 
tions in the UV and X-ray band are significantly 
contributed by the non-thermal radiation from the 
jet, i.e., variations of the synchrotron and SSC 


8 The light travel time from the corona to the disk is generally 
negligible, of the order ~ 10 3 s assuming a scale-height of 
the corona as a few Schwarzschild radii above the disk. 
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components respectively, which are correlated by 
themselves H. 


5.4. Constraining the Black Hole Mass 


Given their relatively small widths of the broad 
lines, the black hole masses in NLSls are believed 
to be systematically lower tha n those i n classi¬ 
cal Seyfert Is and quasars (e.g. iGrune fe Mathur 
2004 feotte et~ahll2004 Xu et al. 2012h estimated 


from the commonly used virial method. This, 
and the resulting high Edding ton ratios (e.g. 
Grupe et all 2 Q 10 UXu et ahll 2012 ) can explain nat¬ 


urally _mostof the observed properties of NLSls 
(see Komossa 20081 for a review). Major prop¬ 
erties of NLSls include their location at one ex- 


space ( ] 

Boroson 

Leighly 

19994 


.. _. . - rapid X -ray variability (e.g. 

19994 lAi et al.112011 ). relatively low 


plitu des of optical/UV variability (e.g. lAi et al 


1201.4 ). strong gaseous out flows as traced by th e 
[OIII]5007 emission (e.g. iKomossa fc Xu 2007 ). 
and the morphology (pseudo-bulges) of their 


1 o S! 

host galaxies (e.g . lOrban de Xivrv et al.l 12011 


Mathur et al.ll2012l) . However, there remains con¬ 


troversy as to whether their black hole masses 
are underestimated, if th eir broad li ne regions ar e 
planar and seen face on ( Jarvis fc McLurel 2006 1. 
Specifically, with argum ents along the same line, 


Calderonc et al. (120131 ) proposed recently that 


radio-loud NLSls have supermassive black holes 
similar to those in blazars (~ 10 8 — 1O 9 M 0 ). 


An independent constraint on the black hole 
masses of AGNs comes from the variability of their 
X-ray emission, which carries the information 
about the dynamics of the hot corona in the close 
vicinity of the black holes. The Suzaku observation 
shows that 1H 0323+342 varied significantly on a 
timescale of ~20 ks (Figure 0J). Short timescales 


9 S hort ly before the submission of our paper, a paper 
by iPaliva et al.l 120141 ) appeared on the astro-ph preprint 
server. Compared to their work, we reached similar conclu¬ 
sions regarding the Swift light curve, but focus on different 
temporal epochs. We note that they reported the Fermi 

y-ray light curve and concluded that the 7 -ray flux of this 
source was nearly stable throughout years until the end of 
2012 , when anomalous variations in all the bands including 
the 7 -ray band were detected. The data we present in this 
work do not cover that period of time. Our result supports 
the idea that 1H 0323+342 likely spent most of the time 
(at least during 2010 October-November) in quiescence, in 
which its UV/optical and X-ray emission is dominated by 
the disk/corona component. 


as such are almost ubiquitous to AGNs w ith rel¬ 
atively small black h ole masses ( Leighlvl 1999a; 
Miniutti et al. 20091 : Ai et al. 201 4 I Yuan et al 


2014I) . This minimum variability timescale sets 


an upper limit on the light-crossing time of the 
X-ray emitting region with a typical size of ~10 
Schwarzschild radii, which corresponds to a black 
hole mass < 4 x 1O 7 M 0 . 

A more quantitative estimation can be achieved 
by using the X-ray normalized excess variance 
of ms , defined as the variance of the light curves 
divided by the mean flux with the contribution 
from the measurement errors subtract ed (e.g. 
Nandra et alJIl997t IVaughan et al.l 12003 ). Inter¬ 


esting enough, a remarkably tight linear relation¬ 
ship between the short-timescale X-ray normal¬ 
ized excess variances and the black hole masses 


was s u ggested in recent studi es bv lO’Neill et al. 
( 2005 ). Zhou et al. ( 201(4) and Ponti et al. (|2012l + 


This provides a useful tool to estimate the black 
hole masses independent of the commonly used 
virial method. To make it directly comparable 
to previous results and also to eliminate any po¬ 
tential contribution from the jet at high ener¬ 
gies, we extract the 2-4 keV light curva 10 ! from 
the Suzaku observation and divide it into four 
segments of 40 ks each. The normalized excess 
var iance and it s unce rtainty is calculated follow¬ 
ing IPonti et al. ( 20121 ). We find an estimation of 
CT rms, 2-4 = 12.3l 3 ;g x 10 -3 , which corresponds to a 
black hole mass of Mbh = 8 . 6+27 x 10 6 M e using 
the relation obtained from the AGN sample with 
black hole masses estimated via the reverberation 


mapping method (jPonti et al.1 12012 ). This value is 
consistent with that estimated from the broad H/3 
line (~ 1.8 x 1O 7 M 0 , Zhou et aLll2007l) . whose un¬ 
certainty is 0.3-0.4 dex. In Figure[l2l we compare 
the obtained of ms and the Mbh estimated from 
the broad line width of H/3 for 1H 0323+342 with 
the log of ms -log Mbh relation. An error bar indi¬ 
cating the typical uncertainty of Mbh (0.3 dex) is 
also plotted. As can be seen, the two methods give 
broadly consistent Mbh values within the uncer¬ 
tainties. We conclude that, if the observed X-ray 
variability is indeed predominantly contributed by 


10 It should be noted that this relation is insensitive to the 
energy band of the X-ray emission within which oj? ms is 
calculated. Extracting a light curve in the 0.2-2 keV band 
gives <r^ ms 0 2_2 = l+lj/iPj 2 x 10 -3 , leading to a slightly 
smaller Mbh. 
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the emission from the corona instead of from the 
jet, as is suggested in this work, 1H 0323+342 
harbors most likely a small black hole with the 
mass of the order of 10 7 Mq. 
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A. Testing the Significance of the Cross-Correlation Function 


The Monte Carlo technique is employed to generate the simulated X-ray light curves. First, we interpolate 
the observed X-ray light curve into bins of 0.1 day and calculate the power spectrum of the interpol ated light 
curve. We generate the red noise light curves according to the method of lTimmer fe Konig ( 1995 1 with the 
obtained power spectrum. These generated light curves with bins of 0.1 day are then sampled in a fashion 
that mimics the real observations. The Poisson noise is added to each point, considering the contribution 
from the observational uncertainty. In this way, the uncertainties which arise from the stochastic property of 
the red noise and the measurement are both included in the simulation. By repeating this process for 2000 
times, we have 2000 simulated light curves which are uncorrelated with the observed one. Then we calculate 
the ICCF between these simulated light curves and the w2 light curve exactly as done for the observed 
X-rays. The values of each ICCF are recorded so that there is a distribution of ICCF values at each r. The 
68 %, 95% and 99% extremes of the distributions are then plotted in Figure [7] as dashed, dash-dotted and 
dotted lines, respectively. 
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Table 1: Results of the UVOT calibration proce- 


d ure as described in Section 12.II 


Band 

^eff 

(A) 

CF 

(ergs -1 cm -2 A -1 ) 

A a 

(mag) 

V 

5427 

2.60 x 10 -ly 

0.65 

b 

4355 

1.46 x 10 -16 

0.86 

u 

3473 

1.64 x 10 -16 

1.04 

w 1 

2611 

4.31 x 10 -16 

1.48 

m2 

2255 

8.30 x 10 -16 

1.79 

w2 

2082 

5.86 x 10 -16 

1.72 
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Table 2: UV and X-ray Variability Amplitudes and 
the Results of Spearman’s Rank Correlation Test 


Band 

Avar 

P 

P 


(%) 



(1) 

(2) 

(3) 

(4) 

w\ 

14.9T0.4 

0.66 

2.8 x 10~ B 

4.8T0.9 

0.77 

0.7% 

m2 

15.7T0.5 

0.75 

1.3 x 10- 7 

7.0T1.1 

0.88 

0 .1% 

w2 

14.3T0.3 

0.75 

< 10 -8 

9.3±0.4 

0.73 

1 x 10" 5 

X-ray 

33.9T0.6 

28.8T0.8 




Note. —The <r m , F var , p and p in the first line of each band are values for the whole dataset of the Swift observations used 
in this paper, and in the second line are values for the observations during 2010 October-November. Column 2 is the fractional 
variability amplitude calculated using equation CD: Column 3 and 4 are the Spearman’s rank correlation coefficient and the 
corresponding p- value. 
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Table 3: Parameters of Four Models Applied to 
the Time-Averaged Spectrum of the Suzaku Data 


Power law+black body 

r 

1.87 ±0.02 

kT [keV] 

0.15 ±0.01 


X+d.o.f. 

349/357 

Double Power law 

Ti 

Norm a r 2 

Norm 

X+d.o.f. 

q 1 q+0.19 
°- ±o -0.17 

2.49 ±0.36 1-58+0® 

i qi+°- 33 
i ' al -0.35 

319/357 

Power law±Partial Covering Absorption 

r 

1V H D logC A c 


X+d.o.f. 

2.40 ±0.04 

6.06±g;lg —0.76 ± 0.63 

0.58 ±0.03 

391/356 

Power law±Comptonization 

r 

1 70+0.04 

1 - /O o.i2 

fcT plasma [keV] r e 

0.28+;®® 13.41+ 


X 2 /d.o.f. 

330/359 

Power law±Reflection 

r 

1.95±0.01 

£r S 

0.92++ 225+® 

< 0.13 

X+d.o.f. 

351/355 


a Flux density of the power-law component at 1 keV in units of 10~ 3 photons keV —1 cm -2 s _1 . 
b Absorption column density TVh of the partial covering absorber in units of 10 22 cm~ 2 . 
c Ionization parameter of ionized partial covering model, defined in I Reeves et alj d2008t) . 
d Covering fraction of the absorber. 
e Optical depth. 

f Abundance of iron relative to solar value. 

^Ionization parameter of reflection model, defined in lRoss Sz Fabianl (I2005T) . 

^Spin parameter of the reflection model. The model cannot constrain this parameter and only gives an upper limit with 90% 
confidence level. 

Note. —The neutral absorption column density is fixed at the Galactic value N^ al = 1.27 x 10 21 cm -2 . 
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Table 4: Parameters of the Absorbed Power Law 
Plus black body Model Applied to the High and 
Low State Suzaku Spectrum 



r 

kT [keV] 

X 2 / d.o.f. 

R a 

High-state 

1 QQ+U.U5 

o.o6 

0.15 ±0.02 

214/232 

11% 

Low-state 

1.80 ±0.06 

0.15 ±0.02 

153/170 

11% 


a The 0.3-10 keV flux ratio of black body to power-law component in the model. 

Note. —The neutral absorption column density is fixed at the Galactic value 7Vp^ al = 1.27 x 10 21 cm -2 . 
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Table 5: The GHz Radio Data of the SED of 
1H 0323+342 Retrieved from NED_ 


Frequency 

(GHz) 

Flux 

(mJy) 

References 

8.4 

277 

1 

5 

357 

2 

4.89 

304 

3 

4.85 

368 

4 

4.85 

362 

5 

4.83 

358 

6 

1.4 

474 

7 

1.4 

614.3 

8 


References. — (1) Healey et al. 2007; (2) Linford et al. 2012; (3) Laurent-Muehleisen et al. 1997; (4) Gregory & Condon 
1991: (5) Becker, White, & Edwards 1991; (6) Griffith et al. 1990; (7) White & Becker 1992; (8) Condon et al. 1998. 
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Table 6: Parameters of the Best-Fit SED Model 


Parameters in the IC/BLR and IC/IR Cases 


Parameter 

(i) 

IC/BLR 

(2) 

IC/IR 

(3) 

7min 

130±15 

28±6 

7b b 

1073±412 

296±90 

7max C 

8000 

12000 

Pl d 

-1.8 

1.8 

P2 e 

5.0 

4.6 

5 { 

2.7±0.6 

13.6±0.9 

B s 

1.9±0.7 

0.34=0.25 

MBH h 

1.8 X 10 * * * * * * 7 

1.8 X 10 7 9 * 

^syn 

(2.0 ± 1.6) X 10 13 

(1.8 ± 1.5) X 10 12 

i j 1 J 

‘•'syn b u 

(1.9 ± 1.7) X 10 44 

(1.1 ± 0.7) x 10 44 

Lbol k 

1.9 x 10 45 

1.9 X 10 45 

-Pd 

7.1 x 10 43 

9.2 x 10 42 

Pe“ 

2.7 X 10 43 

5.4 x 10 43 

Pp 11 

4.6 X 10 43 

9.6 X 10 44 

Pb° 

4.3 X 10 42 

7.7 X 10 43 


“Minimum Lorentz factor of the injected electrons. 
b Break Lorentz factor of the injected electrons. 

“Maximum Lorentz factor of the injected electrons. 

d Slope of the injected electron distribution below 7b- 

“Slope of the injected electron distribution above 7b- 

^Doppler boosting factor. 

9 Magnetic field in units of Gauss. 

h Black hole mass in units of solar mass Mq. 

*Peak frequency of the synchrotron component in units of Hz. 

9 Luminosity at peak frequency of the synchrotron component in units of ergs -1 . 

^Bolometric luminosity Lbol — Ldisk + ^corona in units of ergs -1 . 

1 Radiative power in units of ergs -1 . 
m Power in the bulk motion of electrons of the jet in units of ergs -1 . 

“Power in the bulk motion of protons of the jet in units of ergs -1 . 

“Poynting power in units of ergs -1 . 
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Fig. 1.— Upper six panels: Light curves in the v, b, 
u, w 1, m2, w 2 band, respectively, of Swift /UVOT 
from 2006 to 2013. The fluxes are in units of 
10 -16 ergs -1 cm -2 A -1 . Bottom panel: X-ray 
count rates in the band 0.3-10 keV observed by the 
Swift/X RT. The observations in 2010 October- 
November (in the grey area) are used to inves¬ 
tigate the correlations of the UV and X-ray vari¬ 
ability on short timescales. 
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0.0 0.1 0.2 0.3 0.4 

Count Rate (0.3-10 keV) 


Fig. 2.— The X-ray count rates versus UV 
fluxes for all the Swift observations analyzed in 
this paper. The UV fluxes are in units of 
10 -16 ergs -1 cm -2 A -1 . 
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Days from 2010 Oct 28 


Fig. 3.- Densely sampled X-ray and w2 light 
curves in 2010 October-November. The w2 flux is 
in units of 10” 16 ergs _1 cm -2 A -1 . A dip in both 
the X-ray and w2 band of 1H 0323+342 is present 
around the 10th day. 
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Fig. 4.- The light curves of 1H 0323+342 in dif¬ 
ferent energy bands observed by the Suzaku/X IS, 
from top to bottom: the full XIS band light 
curve (0.2-12keV), the soft band light curve (0.2- 
2keV) and the hard band light curve (2-10keV). 
The bars in the top panel cover the time ranges 
for which the high- and low-state spectra are ex¬ 
tracted. The four 40 ks long segments which are 
used to calculate the normalized excess variance 
are separated by the vertical dashed lines (see Sec¬ 
tion 15.41) . 
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0.1 0.2 0.3 0.4 

Count Rate (0.3-10 keV) 

Fig. 5.- The X-ray count rates versus the w2 
fluxes in units of 10 -16 ergs -1 cm -2 A -1 in 2010 
October-November. The intervals between ad¬ 
jacent observations in the w 2 and X-rays are 1 
3 days. The results of the Spearman rank correla¬ 
tion test can be found in Table [2j together with 
the other two UV bands. 
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x (days) 


Fig. 6. — Cross-correlation functions between 
the X-ray and w2 band during 2010 October- 
November: ICCF (black solid line) and ZDCF 
(dots with error bars). Positive time lags corre¬ 
spond to the X-rays leading. 
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Fig. 7. The ICCF between X-ray and w 2 
band light curves during 2010 October-November 
is plotted as black solid line. Positive time lags 
correspond to X-rays leading. The dashed, dot- 
dashed and dotted lines represent 68%, 95% and 
99% extremes of the distribution of the ICCFs be¬ 
tween simulated X-ray light curves and the ob¬ 
served w2 light curve. 
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1 2 5 

Energy (keV) 


Fig. 8.-- Panel (a) shows the averaged Suzaku FI 
spectrum of 1H 0323+342 (crosses) and the power 
law plus black body model (solid line) fitted to the 
spectrum. Lower panels show the residuals of the 
spectral fits in terms of a with error bars of size 
one. Panel (b): Residuals of single power law fit¬ 
ted to the spectrum. The positive residuals are re¬ 
vealed below IkeV. Panel (c): Residuals of power 
law plus black body model fitted to the spectrum. 
Panel (d): Residuals of the ionized partial covering 
absorption model fitted to the spectrum. Panel 
(e): Residuals of reflection model. The data be¬ 
tween 1.7-2.3keV are ignored on account of pos¬ 
sible calibration uncertainties (see Section [272j) . 
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Fig. 9.- - 1H 0323+342 FI spectra and their best- 
fit power law plus black body models in high-state 
(red) and low-state (green), respectively, during 
the Suzaku observation. See Figure 0] for the 
time intervals of high- and low-state. The spec¬ 
tral shape does not change, though the flux var¬ 
ied. See the electronic edition of the Journal for a 
color version of this figure. 
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Days from 2010 Oct 28 


Fig. 10. The variations of the photon indices 
T of the Swift/X RT observations during 2010 
October-November. The spectra are fitted by a 
single power law with the absorption column den¬ 
sity fixed at the Galactic value. 
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Fig. 11.— The rest-frame broad band SED of 
1H 0323+342. The red circles represent the in¬ 
frared data from the A11WISE Source Catalog. 
The black arrows indicate the 90% upper limits 
from the Suzaku/HXD and the red arrows indi¬ 
cate the 2-er upper limits from the Fermi/ LAT. 
The UV/optical data are corrected for reddening 
by Galactic dust and the X-ray data are corrected 
for neutral interstellar absorption. The radio data 
are not involved in the SED fitting. The green, 
orange and black dashed lines represent the syn¬ 
chrotron, SSC and EC components, respectively. 
The dotted lines refer to the multi-temperature 
black body of a standard disk. The blue solid 
curve is the sum of all contributions from the best- 
fit model. The black hole mass employed here is 
1.8 x 10 7 Mq. The grey solid lines represent the 
SED model fitting results after the hard X-rays 
were varied by a factor of two. See the electronic 
edition of the Journal for a color version of this 
figure. 
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Fig. 12.— The log MaH-log af ms relation cali- 
brated b y t he re verberation-mapped sample in 
Pont! et ah ( 2012h (solid line). The filled circles 


indicate sources in that sample and the open cir¬ 
cle is 1H 0323+342. The errors on the normalized 
excess variances are at 1 -a confidence level. A typ¬ 
ical uncertainty (0.3 dex) of the black hole mass 
measured by the virial method is also plotted for 
1H 0 323+342. The sou rces with only upper lim¬ 
its in lPonti ct al. ( 20121) are not shown here. The 
black hole masses of sources in the reverberation- 
mapped sample were determined using the rever¬ 
beration mapping method. 
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